Complete light absorption in graphene-metamaterial corrugated structures 
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We show that surface-plasmon polaritons excited in negative permittivity metamaterials having 
shallow periodic surface corrugation profiles can be explored to push the absorption of single and 
continuous sheets of graphene up to 100%. In the relaxation regime, the position of the plasmonic 
resonances of the hybrid system is determined by the plasma frequency of the metamaterial, allowing 
the frequency range for enhanced absorption to be set without the need of engineering graphene. 

PACS numbers: 78.67.Wj, 73.20.Mf, 42.25.Bs, 78.20.Ci 



I. INTRODUCTION 

The presence of a boundary in a metal allows for sur- 
face collective charge density oscillations in addition to 
volume plasmon modes. 1 When these surface modes 
couple to light, hybrid photon-plasma excitations emerge, 
known as surface-plasmon polaritons (SPPs).[2H5] 

A notable feature of a SPP is its small wavelength, 
when gauged against the wavelength of a transverse pho- 
ton with the same frequency. This results in high de- 
grees of electromagnetic energy confinement, endowing 
SPPs with high sensitiveness to surface conditions. [H 
[B] For this reason, SPPs are responsible for surface- 
enhanced Raman scattering, [7] enabling single molecule 
detection, [5] and having many technological applica- 
tions (e.g., chemical and biological sensors, [9] and nano- 
resolution imaging|10j). Moreover, the compact nature of 
SPPs offers the prospect of combining large bandwidth 
and integration, at subwavelength scale, eventually allow- 
ing nanoscale photonic circuit implementation. [101 111) 

The scope of surface plasmon-related phenomena 
has been increasing swiftly in the past decade, span- 
ning metallic nanowires,|T2i bandgap nanostructures. [T3] 
metallic nanoparticles,[Tl] and other nano-engineered 
materials. [15HT7] With the advent of truly two- 
dimensional crystals — graphene being the first of 
the kind[IHl EH — a new playground for plasmonics 
has emerged. [20] The unique electronic properties of 
graphene, characterized by massless and chiral low- 
energy excitations, [5T] and unconventional transport 
properties, [22 originate the most distinct electromag- 
netic confinement behavior: SPPs with propagation 
lengths exceeding those of conventional metal-dielectric 
interfaces, (23j and guided transverse-electric modes. [24l 
125] just to mention a few. 

The observation of prominent plasmonic absorption 
peaks in graphene microarrays[26 has triggered a new 
research line, [27J HH] m which plasmonic excitations are 
explored to overcome the major obstacle in graphene- 
based optoelectronics: the small light absorption in one- 
atom thick graphene samples. 29 Inspired by these works 



and well-established results for metallic gratings, [5] we 
propose a hybrid graphene-metamaterial system, where 
a single and continuous graphene sheet is seen to absorb 
all the light impinging on it. The frequency range for 
enhanced absorption is determined by the metamaterial 
alone (see later), which allows subsequent control over 
the absorbed plasmonic waves to be performed in the 
graphene sheet (e.g., via chemical doping[30j. 




Figure 1: Schematic picture showing transverse magnetic light 
impinging on graphene placed on a periodically structured 
conducting substrate. The light is converted into surface plas- 
monic waves and subsequently absorbed in graphene. 

The system we have in mind is depicted in Fig. [TJ inci- 
dent light with angular frequency uj strikes the interface 
S of a periodically structured conducting system coated 
with graphene. The light is transverse magnetic polar- 
ized, H in (r,t) = H in e i( - k - r - ut )e y , where k = (k x ,0,k z ), 
and the parametrization of the boundary reads X = 
{(x,y,z)\z — ((x)}, with £(x) having the symmetry 
property, C( x + a) = ((%)■ I n order to achieve electro- 
magnetic energy confinement and SPP-related phenom- 
ena at the interface S, the substrate is assumed to have 
e(u) < in the frequency domain of interest. [5J In a prac- 
tical implementation, the latter would be an engineered 
metamaterial. [3T] in the relaxation regime uj p 3> u> 3> 7 
(with uj p being the effective plasma frequency and 7 the 
relaxation rate). 

The spatial periodicity of the present configuration is 
instrumental in order to excite SPPs via coupling to the 
radiation field: conservation of momentum requires that 
the projection of the wavevector of the incident photons 
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onto the plane parallel to £ matches the wavevector q 
of SPPs belonging to the hybrid structure. In flat sur- 
faces the matching condition cannot be achieved since, 
in that case, the SPP dispersion relation lies to the right 
of the light line (i.e., q > ui/c > k x ). The role of the 
periodic corrugated profile is thus to open Bragg scatter- 
ing channels providing photons with the necessary wave 
vector (according to q — k x ±2Trn/ 'a, for some integer n), 
thus allowing the excitation of SPPs to occur via direct 
illumination. 



II. LIGHT SCATTERING: EXACT 
INTEGRAL-EQUATION FORMULATION 

The investigation of the optical properties of the hy- 
brid graphene-metamaterial structure (Fig. [T]) is carried 
out via exact numerical calculations of the reflected and 
transmitted fields. The total magnetic field is given 
by H v (r,t) = H±(r) e - luJt [here r = (x,z), and ± = 
sign (z — £(x))], where the amplitudes H±(r) are solu- 
tions of the wave equation, [V 2 + e±(w)w 2 /c 2 ] H±(r) = 
0, with V 2 = d 2 +d 2 , and where e+(w) = 1 (vacuum) and 
e_(u;) = 1 — ujp/[uj(u! + ij)] (substrate). The boundary 
conditions at the interface £ read as 32 
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where the notation means that the fields (and deriva- 
tives) are to be evaluated taking values of z approaching 
£ from above (C + ) or below (£~) , d n = n ■ V denotes 
the derivative along the unit vector normal to the surface 
n, directed from the vacuum into the substrate, and J s 
is the longitudinal current induced by graphene. 

The surface current relates to the electromagnetic 
fields according to the constitutive relation (Ohm's law), 
Je = J(x) = a(x,u>)E(x,C) ■ t, where a(x, uj) is the 
optical conductivity of the corrugated graphene sheet, 
which in general may depend on the x coordinate, and 
E(x, C) • t = iL(x) / (soui) is the tangential component of 
the electric field evaluated at r = (x, ({x)). Here, t de- 
notes the unit vector tangent to £ and Eq is the vacuum 
permittivity. In what follows, we approximate the op- 



tical conductivity by its bare value cr(w) in the absence 
of corrugation. This is justified since the inclusion of a 
spatial modulation in cr(uj) does not significantly alter 
the optical properties to be discussed throughout (see 
Sec.|V]). We take c(w) as obtained by the random-phase 
approximation, 22] with a relaxation rate T = 2.5 meV, 
in consistency with typical values. [33H3S] 

In order to solve for the full electromagnetic field we 
make use of an exact integral-equation method devel- 
oped by Toigo et al. [35] This method exploits the in- 
tegral form of Maxwell equations in order to express 
H±(r) as an integral over the functions H(x) and L(x) = 
L(x)y/l + [('{x)} 2 . Making use of the Bloch prop- 
erty, H(L) = En H ( L )ne iqnX , with q n = k x + 2im/a, 
Toigo and co-workers reduce the well-known integral 
equations [37] to a set of linear algebraic equations for 
the Fourier coefficients of H±(r). Here, we generalize 
their method as to include the effect of a metallic sheet 
with conductivity cr(w) deposited on the corrugated sur- 
face £. We just state the basic results: above (below) 
the selvedge region z > max C, (z < min £) , the re- 
flected (s = +) [transmitted (s = — )] wave is given by 
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responds to the specular wave, and the remaining chan- 
nels correspond to propagating modes (Bragg beams) or 
modes confined to the surface £ (evanescent waves), de- 
pending on whether 
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is pure imaginary or real, respectively. The amplitudes 
T-L s m relate to the Fourier coefficients {H n , L n } of the func- 
tions defined in Eqs. ([T]) and ^ according to 
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where T^f(p) = a" 1 /°^ 2 d xe -2*W<x e -a/3K„(g m )c(z) 

(a,/3 = ±1) encodes the effect of the spatial profile of 
the interface, and 
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provides an infinite set of linear equations for the un- tions is solved by retaining a finite number of coeffi- 
known coefficients {H n , L n }.[3H] This system of equa- cients, n = —N+ 1...N. The integer N must be chosen 
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sufficiently large so that convergence is obtained. The 
actual convergence properties of the present method de- 
pend crucially on the form of C(x), and, in general, for 
corrugations with very large amplitudes and/or possess- 
ing discontinuities, convergence is not guaranteed. pM |4"0] 
Here, we focus on periodic structures with a sinusoidal 
profile, C( x ) — Co cos (2"Kx/a), for which the functions 
T^f (p) have a well-known form and convergence is fast, 
i.e., only a few Fourier coefficients are required for good 
accuracy. [H] 



III. MAIN RESULTS 

We focus our discussion on normal incidence (9 = 0), 
with the important remark that the reported features 
can be observed at 9 > 0. In that case, incident pho- 
tons have k x = 0, and hence the wavevectors of SPPs 
belonging to the hybrid structure coincide with the re- 
ciprocal lattice wavevectors, q n = G n = 2n/a. In the 
relaxation regime and in the ideal scenario of no losses 
occurring inside the metamaterial (except when stated 
otherwise, we assume 7 = 0), the fields cannot propagate 
into the substrate, and hence the incident light is totally 
reflected at the interface (in the absence of graphene). 
Yet, signatures of coupling of the incident light to SPPs 
are observed in the zeroth order (specular) reflectivity of 
metallic gratings. [6l |42l |43] The signatures arise because 
surface plasmonic waves excited at vacuum-metal struc- 
tured interfaces via direct illumination are radiative: the 
excited SPPs reradiate into the vacuum via open Bragg 
channels (w > c|G„|, with |n| > 0). The resulting gen- 
eral decrease in the specular reflectivity is accompanied 
by pronounced features (e.g., dips and bumps) appearing 
at the SPPs resonances. 

Coating the metamaterial interface with graphene in- 
troduces new channels for SPP dissipation, as encoded 
in the optical conductivity ct(cj), due to Drude diffu- 
sive scattering and interband transitions. [52] Based on 
general arguments, [44] we expect that if the SPP decay 
rate, r g , introduced by graphene, matches the broaden- 
ing due to radiative decay, r ra d, the interference between 
the specular channel and the high-order Bragg beams 
suppresses the reflectivity, hence allowing the SPPs to 
be fully relaxed within the graphene sheet. The latter is 
borne out in Fig. [2] (right top panel) , which shows the ab- 
sorption near the lowest SPP resonance for a sinusoidal 
corrugation with £o = 85 nm and a — 8 /im, and at fixed 
Fermi energy, Ep = 50 meV. We note that by varying 
the electronic density of the graphene sheet, it is pos- 
sible to obtain complete absorption for Co/a as low as 
10~ 3 (see Fig. [3|. We thus predict that small corruga- 
tion amplitudes are needed for obtaining complete light 
absorption. 
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Figure 2: The SPP spectrum for a flat-surface interface E 
with huj p = 0.5 eV and 7 = (solid line in the left panel). 
Its intersection with the displaced incident photon line (k x — > 
k x + G n ) gives an estimate of the frequency of SPPs in the 
nth band [e.g., Ml(Gi) ~ 0.147 eV]. Absorption at normal 
incidence near the first SPP resonance (right top panel), and 
the scattering function S(u) for a wider frequency interval 
(right bottom panel). Dashed lines in the right panels refer 
to the bare interface (in the absence of graphene). 
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Figure 3: Dependence of the absorption at resonance on the 
corrugation height (left panel) for different graphene Fermi 
energies, Ef, and relative absorption for a nonzero metama- 
terial relaxation's rate 7 (right panel) with 7/cjp = {10 -4 , 5 • 
10 -4 , 10" 3 } (magnitude increasing as indicated by the arrow) . 
(Left panel) The fast change in A(Gi) as Ef varies from 0.07 
to 0.08 eV results from a suppression of interband transitions 
in graphene as 2Ef > h£l(Gi). (Right panel) Appearnce of 
two peaks for the highest values of 7 witnesses a considerable 
shift of the resonant position of A(Gi) relative to ^lbare(Gi). 
Other parameters as in Fig. [2] 



IV. DISCUSSION 

In order to understand the reported SPP-assisted en- 
hanced absorption, we first characterize the spectrum 
of SPPs in the hybrid graphene-metamaterial system. 
The effect of a periodic spatial profile can be qualita- 
tively appreciated by analogy with one-dimensional elec- 
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trons subjected to a weak periodic potential: a length 
scale is introduced (a), and the SPP spectrum acquires a 
multiband structure with period 2ir/a, with gaps open- 
ing at the Brillouin zone boundaries ±7m/a.[4"5] In- 
deed, SPP modes can be excited by photons via cou- 
pling through lattice reciprocal vectors, G n = 2-nn/a, 
with \n\ > 0, whenever the frequency of the incident 
photon matches the frequency of a SPP in a given 
band, that is, uj = ft(q n ). For shallow spatial pro- 
files, the SPPs frequencies can be estimated by folding 
the flat-surface spectrum in the absence of graphene. [4"5] 
k = (n/c)^e(n)/[e(Q,) + 1] [with e(O) < -1], into the 
Brillouin zone and locating the intersection with the inci- 
dent photon line [2j |3j 0143] (or equivalently, by the proce- 
dure illustrated in the left panel of Fig. [2]). For a rigorous 
determination of the SPP spectrum for arbitrary corru- 
gation profiles, we examine the poles of the scattering 
matrix S, defined by (U + ,U^) T = S(H in ,0) T , where 
"H* = ('H ± N+1 ,...,'H%) and H in are 2iV-dimensional 
vectors containing the incoming (outgoing) modes in 
the vacuum (substrate), and the input field amplitudes, 
H™ — Hi n 5 n! o, respectively. The SPPs frequencies for 
a given wavevector k x can be obtained by localizing the 
peaks of the function S(w) = m I I • @Z] The scat- 
tering function for normal incidence (k x = 0) is given 
in the right bottom panel of Fig. [2] Clearly, the SPPs 
frequencies are not significantly altered by the presence 
of graphene (as seen by the agreement of the peaks posi- 
tions in both the solid and dashed lines |46j). Moreover, 
these frequencies agree fairly well with the flat-spectrum 
estimate, shown in the left panel, based on perturba- 
tion theory arguments. The effect of graphene is there- 
fore to change the spectral weight of scattering function 
S(lu) around the vacuum-metal resonances, lo ~ f2(G n ), 
without affecting too much their position (note that in 
Fig. [2]only the first two resonances, \n\ = 1,2, are shown). 
The little sensitiveness of the resonant frequency of SPPs 
to the graphene sheet characteristics turns out to be of 
significance for experimental realizations. We can imag- 
ine that by changing e(u>), for example, by considering a 
metamaterial with an appropriate plasma frequency, the 
resonant frequencies £l(G n ) can be controlled within a 
broad range of values, thus tailoring the hybrid structure 
for a desired spectral region. Another possibility is to 
probe the Fermi energy of graphene by optical measure- 
ments, either by monitoring changes in the reflectivity or 
by measuring small changes in the SPPs frequencies. 

The dependence of the absorption line-shape on the 
corrugation amplitude shown in Fig. [3] (left panel) is 
similar to that found in metallic nanopatterncd gratings 
with considerable ohmic losses. [131 HI] Initially, as Co in- 
creases, so does the coupling of SPPs to the radiation 
continuum, hence originating a monotonous increase in 
the absorption at resonance, A(G\). At some point, the 
radiative decay balances the SPP damping in the vicinity 
of graphene, r ra d = r g , and the reflectivity is completely 
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Figure 4: Contour plot for the magnetic field modulus in the 
hybrid graphene-metal corrugated interface (top) and for the 
same structure in the absence of graphene (bottom). The field 
modulus is periodic in x and thus only the restricted interval 
\x\ < a/2 is shown. The corrugation height is £o = 85 nm and 
the impinging light has fkv — 0.1474 eV. The bar indicates 
the scale of the oscillation resulting from interference of the 
incident field and the reflected field, 2Ai n = 4-kc/lj. 

suppressed (absorption is maximum) . [441 |4"7] Increasing 
Co beyond such value overcouples the SPP to the radi- 
ation field, and the specular wave is no longer canceled 
out by the reradiated field: A{G\) decreases. Eventually, 
r ra d r g , and absorption becomes negligible, with all 
the energy being carried by the reflected wave. Evidently, 
the value of Co leading to complete light absorption de- 
pends on the Fermi energy of graphene, since the latter 
determines the effective value of r g . For instance, for 
Ep = 0.1 eV, the maximum absorption occurs for very 
small Co (see Fig. |3|, and given that, at lowest order, 
Trad ~ Co ( see J e -S-i Ref. HSJ, we can infer that the SPP 
damping in graphene T g is particularly small in such a 
case. 

The calculation of the field intensity distribution in the 
vicinity of the surface £ further confirms the role played 
by the SPP decay channel introduced by graphene: In 
Fig. [4] we plot the magnetic field intensity in the vac- 
uum for incident light in resonance with SPPs in the 
first band, \n\ = 1. When graphene is absent (bot- 
tom panel), the large field enhancement effect typical of 
SPPs in corrugated interfaces is clearly observed, [43] 08] 
with \H y \ ss 60 1 i?i n I in the near field. Here, the reflec- 
tivity is total as inferred by the sharp field interference 
pattern observed for z 3> Co- When graphene is intro- 
duced (top panel), the far-field interference pattern dis- 
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appears, demonstrating the suppression of the reflected 
wave, while the near field displays more modest field en- 
hancements, due to SPP damping in the graphene sheet, 
in consistency with arguments given previously. 

V. FINAL REMARKS 

We finish our discussion with a brief analysis of how 
losses in the metallic substrate and modulation of the 
graphene conductivity may affect the optical properties 
of the graphene-metamaterial hybrid system. 

Losses in the metamaterial. When the dieletric func- 
tion of the substrate e(u) is assumed to be a real function, 
absorption in the graphene sheet reaches 1 00% efficiency 
when adequate choices of the corrugation height are made 
(refer to the left panel in Fig. [3p _ On the other hand, in 
any realistic scenario, a trade-off between absorption in 
the graphene sheet and damping within the engineered 
substrate is expected to happen. In order to quantify the 
impact of a non-zero damping parameter 7 (and hence a 
complex dieletric function) on the absorption efficiency 
of graphene paired up with the metallic substrate, we 
compute the change in the absorption, 5 A = A — A^ are . 
Here, Ab are denotes the absorbance of the bare meta- 
material interface (i.e., without graphene). (Note that 
subtracting Ab arc to the absorbance of the hybrid system 
is necessary in order to correctly infer the amount of en- 
ergy absorbed by the graphene sheet.) The right of Fig. [3] 
shows 8 A for a specific configuration (Ep = 0.05 eV and 
Co = 85 nm) and a few values of the damping rate 7. In 
these cases, and also more generally, we have found that 
damping frequencies as low as 7 < 10~ 3 w p are required 
so that the graphene sheet is able to absorb considerable 
amounts of light (SA > 0.1). With this respect, we recall 
that in the infrared range, ratios 7/wp of the order of 10 -3 
are achievable in metals as copper, [33] and hence the re- 
striction on the substrate losses should not impose any 
fundamental limitation on the plasmonic performance of 
the hybrid system. 

Modulation of the conductivity. When graphene ad- 
heres to a corrugated surface its local optical response be- 
comes sensitive to local perturbations (e.g., via position- 
dependent strain). However, in the present case, the SPP 
dispersion relation of the hybrid system is practically un- 
affected by graphene's intrinsic conductivity. |46| and thus 
the impact of modulations in o~{uj) in the optical prop- 
erties of the hybrid system is expected to be weak. In 
order to assess this assumption, we have considered a 
conductivity profile a(x,uj), with the periodicity of the 
corrugated surface and magnitude variations of the or- 
der of 100%, namely, a(x,uj) = cr(w)[l + cos(27nc/a)], 
and no qualitative changes relative to the unmodulated 
case have been found; for instance, the peak of the ab- 
sorbance curve A = A(w) in Fig. [2] is seen to decrease by 
only about 0.1% and shifted by even a lesser amount. 



VI. CONCLUSION 

We have shown that complete light absorption occurs 
in continuous sheets of graphene placed on top of sim- 
ple plasmonic nano-structures. The physical mechanism 
on the basis of the two orders of magnitude boost in 
the graphene's capability to absorb light demonstrated 
here lies in the efficient excitation of SPPs in the pe- 
riodically corrugated interface of the hybrid graphene- 
metamaterial system. Absorbances far exceeding the 
graphene's bare value (of about 2%) [29] are still achiev- 
able when intrinsic losses in the metamaterial are taken 
into account, provided that the damping rate is much 
smaller than relevant frequency scales, namely, the fre- 
quency of light and the plasma frequency of the metama- 
terial. In this so-called relaxation regime, 7 C w < w p , 
the plasmonic resonances of the hybrid system are de- 
termined by the plasma frequency of the metamaterial, 
allowing the frequency range for enhanced absorption to 
be set without the need of engineering graphene. Fur- 
thermore, the plasmonic-enhanced absorbance effect is 
shown to occur in a wide range of electronic densities of 
the graphene sheet (either low- or high- doped samples). 
The properties of the hybrid graphene-metamaterial sys- 
tem described in this work open interesting possibilities, 
such as the capture of light and further guidance and 
control of the excited plasmonic waves (e.g., by design- 
ing SPP propagation paths by inhomogeneous doping 
of the graphene sheet |50j ) . making this system a poten- 
tial candidate for plasmonics applications, such as light- 
harvesting and photonic circuit implementation. 
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